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Abstract 

We compute the next-to-leading, 0(ma 5 ), contribution to the hyperfine splitting in 
positronium within the framework of NRQED. When applied to the ground state, our cal¬ 
culation reproduces known results, providing a further test of NRQED techniques. Besides 
providing a very simple method of calculation of the standard result, we also obtain new 
expressions for excited states of positronium with negligible additional effort. Our cal¬ 
culation requires the complete next-to-leading matching of the lowest-dimension NRQED 
four-fermi couplings, which we publish here for the first time. 


1. Introduction and Summary 


The detailed comparisons between the heroic precision calculations using Quantum 
Electrodynamics (QED) and equally heroic precision measurements represent one of the 
pinnacles of twentieth-century physics. On the theoretical side, comparisons involving 
the energy levels of bound states, such as hydrogen or positronium, pose a particular 
challenge. This is because of the necessity of incorporating the small radiative corrections 
of relativistic QED into the nonperturbative treatment required for the bound state itself. 
Important conceptual progress in handling bound states in QED was made several years ago 
[1] with the development of nonrelativistic QED (NRQED), which consists of an application 
of effective-field-theory ideas to atomic physics applications of QED. 

NRQED starts with the recognition that much of the complication of QED bound-state 
calculations arises because the fine-structure constant, 1 a = e 2 /47r, enters into observables 
in two conceptually different ways. First, a enters as the small parameter which controls 
the higher-order QED radiative corrections. Second, a enters from the appearance in these 
calculations of three separate scales: the electron mass, m, the bound-state momentum, 
mu, and the bound-state energy, mv 2 , once it is recognized that v = 0(a) in the bound 
state. While much of the higher-order QED radiative corrections involve scales at, or 
above, m, where relativistic kinematics is important, the complications associated with 
handling the bound states all arise at the lower two scales, mv and mv 2 . 

NRQED takes advantage of this hierarchy of scales to efficiently separate the radiative 
corrections from the bound-state physics. First one accurately integrates out all physics 
associated with momenta p > 0(m ), obtaining an effective theory of nonrelativistic parti¬ 
cles whose interactions are organized according to their suppression by powers of the two 
independent small parameters, a and v. This effective theory is then used to systemat¬ 
ically compute bound-state properties, at which point v becomes of order a. Keeping v 
and a independent until this last step makes the bookkeeping more straightforward. The 
additional bonus is that bound state calculations are much easier to do within the effective 

1 We use throughout units for which K=c— 1. 


2 



theory because its nonrelativistic framework permits the direct application of well-tested 
techniques based on Schrodinger’s equation. 

Any effective field theory relies on the existence of a hierarchy of scales, say Mi <C M 2 , 
in a physical problem. The heart of the effective theory’s utility lies in its powercount¬ 
ing rules, which identify how to systematically isolate the complete contributions to any 
observable to any fixed order in the small ratio, M 1 /M 2 . NRQED is no exception in this 
regard, with the powercounting rules identifying the suppression of observables in powers of 
v ~ a. It is the recent development of NRQED power-counting rules [2], which now makes 
it possible to directly identify the 0(a n ) contributions to any atomic-physics observable. 
These rules have recently been demonstrated in practice in an illustrative calculation of 
the Lamb shift [3]. 

Although contributions order by order in a can also be obtained by other methods, 
the virtue of NRQED lies in its simplicity, since this potentially brings more complicated 
calculations into the domain of the feasible. Besides permitting the very simple extension of 
low-order results to bound states involving particles of other spins [4], NRQED’s simplicity 
has very recently been exploited to give the first-ever analytical calculation of the last 
previously uncomputed 0(ma 6 ) contribution to positronium hyperhne splitting [5]. 

In the present paper we continue this process in other directions. Since the power of 
effective field theory methods are best seen once one goes beyond leading order in small 
parameters, we present here an NRQED calculation of the complete contributions to the 
positronium hyperhne structure at next-to-leading order, ()(rno~ > ). The result we obtain 
for the ground state hyperhne splitting agrees with previous results [6]. However our results 
are easily extended to the hyperhne splitting of excited states having arbitrary quantum 
numbers n and £. These agree with previous results [7] for arbitrary n but £ = 0. but to 
our knowledge ours is the first calculation which is applicable to general n and l. 

Our goal in presenting this calculation is twofold. One of our aims is to make NRQED 
accessible to a wider community. Detailed calculations which can be compared to known 
results partly act as a vehicle for explaining how NRQED works. O 11 the other hand, 
some of our results are new and interesting in their own right. Our new results come in 
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two forms. First, an intermediate step in our calculation is the matching from QED to 
NRQED at next-to-leading order in a. Only part of this matching is already available 
in the literature, and the remainder, which we publish here for the first time, will have 
applications to many other higher-precision NRQED calculations. In addition, our results 
for the hyperfine structure of the excited states of positronium are also new. 

We organize our presentation as follows. In §2 we briefly review the NRQED la- 
grangian, including a discussion of the matching which is required to obtain some of the 
NRQED couplings to next-to-leading order in a. We first summarize the next-to-leading 
matchings which are already given in the literature, followed by a calculation of a final class 
of matchings which have not been previously published, but which are required to obtain 
the hyperfine structure to the order we require. §3 gives a quick summary of NRQED 
powercounting, as obtained in ref. [2]. These powercounting rules are then used to system¬ 
atically identify all possible contributions to the 0(ma 5 ) hyperfine splitting. The main 
result of this section is that the complete 0(ma 5 ) contribution is obtained in NRQED 
using precisely the same graphs as for the 0(ma 4 ) contribution, but with coupling con¬ 
stants which are matched to relativistic QED to higher order in a. These results are finally 
brought together in §4, where we present our results for the hyperfine splitting. We do 
this for both the ground state and for the excited states of positronium. 

2. NRQED 

We start with the lagrangian density of NRQED as applied to nonrelativistic electrons 
and positrons, which is obtained from full QED by integrating out all virtual physics at 
scales greater than the electron mass, A > m. 

2.1) The Lagrangian 

The fields representing the low energy degrees of freedom are ip, \ and vl /( , which 
respectively represent the nonrelativistic electron, the nonrelativistic positron and photons 
with energy less than m. The lagrangian is: C = Tphoton + ^ 2 -Fermi + ^ 4 -Fermi + • • •, with: 
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Here Dt = dt — iqeAo and D = V — iqe A, where g = — 1 for electrons and g = +1 for 
positrons, a represents the usual Pauli spin matrices. The Feynman rules corresponding 
to this lagrangian are listed in Figure (1). 


2.2) Matching at Leading-Order 


The various coefficients, c t . appearing in this expression are calculable functions of 
a and m which are obtained by integrating out scales larger than m using QED. They 
may be conveniently determined by computing scattering processes for free electrons and 
positrons using both full QED and the NRQED lagrangian, Eqs. (1). Equating the results 
to a fixed order in a and v completely determines the constants c, to this order. 

Performing this operation at tree level in QED gives the lowest-order results in a [1]: 
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We omit the four fermion terms proportional to cr and cs which are listed in ref. [2], because they are 
redundant in the sense that they may be expressed in terms of those we display. 
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Figure (2) 

The matching which determines c^\ 

To see how these arise consider, for example, the contribution to the effective four- 
fermi couplings C4 and C5. These four-fermi operators reproduce in NRQED the effects 
of the tree-level s-channel annihilation graph of QED, given in Figure (2). This must 
be reproduced by an effective interaction because the photon which is exchanged must 
necessarily involve four-momenta of order m, and so cannot appear in the effective theory. 

In performing the matching, we evaluate all diagrams at threshhold — i.e. with van¬ 
ishing external three momentum in the center of mass frame of the electron-positron pair. 
We use this choice of external momentum purely for convenience, and the matching could 
equally well be performed with non-zero (but nonrelativistic) three-momenta for the elec¬ 
tron and positron, giving precisely the same NRQED coefficients. Notice also that care 
must be taken that the scattering amplitude is normalized in the same way when compar¬ 
ing the QED and NRQED results, since the covariant normalization of relativistic spinors 
involves additional factors of \[2E compared to nonrelativistic applications. 

The contributions to the coefficients C 4 and C 5 are distinguished from one another 
by separately evaluating the QED amplitudes for the electron-positron spins combined 
into a spin triplet (S = 1 C4) or a spin singlet (S — 0 <-»• C5). It follows that must 

vanish since the spin-singlet combination of the electron and positron is forbidden by charge 
conjugation invariance from receiving a contribution from the s-channel annihilation graph 
of Fig. (2). 

Although the f-channel photon-exchange graph is also involved in electron-positron 
scattering at tree level in QED, it does not contribute to either c^ ) or c^K Indeed, for 
nonrelativistic electrons and photons the energy of the exchanged photon is much below the 
electron mass and so this scattering is described by the same f-channel graph in NRQED. 
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As a result tree-level t-channel photon exchange does not contribute to any of the NRQED 
four-fermi operators in the matching process. 

It bears emphasis that the matching calculation does not involve any bound states at 
all, since it is done using only the scattering of free electrons and positrons. Matching is 
also the only stage of the calculation which involves QED diagrams. Once the couplings 
of the NRQED lagrangian are obtained in this way, only they are used in the subsequent 
bound state part of the calculation. This separation of the matching from the bound state 
calculations lies at the heart of NRQED’s simplicity. For example, this separation permits 
the use of different gauges in the relativistic and nonrelativistic parts of the calculation, 
since the gauge choice when using NRQED is independent of the gauge used in QED, so 
long as one computes only gauge-invariant quantities. This permits the convenience of 
using a covariant gauge, like Feynman gauge, in the QED part of the calculation, while 
keeping Coulomb gauge for NRQED calculations. 


2.3) Matching at Next-to-Leading Order 


More accurate calculations, such as those of interest in this paper, require the next- 
to-leading corrections to the coefficients of Eqs. (2). Some of these corrections are already 
given in the literature, and some we compute here for the first time. Since our new 
calculations are for the coefficients of the four-fermi coupling constants, C4 and C5, we 
discuss the corrections to these quantities in some detail, and give only a brief summary of 
the higher-order corrections to both £ p hoton and £2-Fermi- We write c t = + ... 

with cf as given in the previous section, and now concentrate on computing the next 
corrections, 


Among the simplest higher-order corrections to the NRQED lagrangian, Eq. (1), are 
the contributions to £ p hoton- To lowest order these are produced by vacuum polarization, 
which gives: 
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Similarly, QED one-loop vertex corrections modify the couplings in £ 2 -Fermi; to give 
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[8], [9]: 
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Here A is the ultraviolet cutoff used to regulate NRQED loop graphs which arise when 
matching. In any calculation of physical properties, such divergences cancel amongst 
themselves, or against explicit cutoff dependence which arises from divergences in NRQED 
loop graphs. 

We now turn to the next-to-leading corrections to -Fermi ■ These corrections may 
be divided into the following three classes according to the topology of the one-loop QED 
graphs which are involved: 

• One-Photon Annihilation: These corrections consist of QED graphs which describe one- 
loop corrections to the tree-level process of the s-channel exchange of a single virtual 
photon. As before f-channel exchange of a single virtual photon does not contribute cor¬ 
rections to T4— Fermi • 

• Two-Photon Annihilation: These consist of the QED ‘box’ graphs which describe the 
s-channel exchange of two virtual photons. 

• t-Channel Two-Photon Exchange: The final class consists of QED box graphs which 
describe the f-channel exchange of two virtual photons. Although f-channel one-photon 
exchange does not contribute to NRQED four-fermion interactions, f-channel two-photon 
exchange does contribute because there is a region of phase space for which the loop mo¬ 
mentum is larger than the electron mass, and so which does not appear in the corresponding 
two-photon exchange graphs of NRQED. The corresponding physics is therefore put back 
into the effective field theory through the NRQED four-fermi interactions. 

We now describe, in turn, the matching due to each of these classes of graphs. While 
the contributions to c\ l > and due to the first two of these may be found in the literature 
[8], [9], those due to the third class we present here for the first time. 
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Figure (3) 

The graphs which give the one-photon annihilation contribution to c^\ 

Matching due to one-loop corrections to single-photon annihilation are described by 
the graphs of Figure (3). The contributions to the scattering amplitude of electrons and 
positrons at threshhold from the QED graphs on the left-hand-side of the equality in this 
figure must be equated to the contributions of the NRQED graphs on the right-hand-side. 
Separately evaluating the QED graphs for a spin-singlet and spin-triplet e + e“ state, and 
solving the resulting equalities for and then gives [ 8 ], [9]: 

(I -7 ann) = 4 ^( 1-7 ann) = 0 . ( 6 ) 

Just as for tree-level matching, the contribution to the spin-singlet coupling, C 5 , vanishes 
by virtue of charge-conjugation invariance. 

Although all possible one-loop s-channel single-photon exchange QED graphs appear 
on the left-hand-side of Fig. (3), the same is not true of the NRQED graphs appearing 
on the right-hand-side. A 11 important issue is therefore how to determine which NRQED 
graphs must be included to any given order in the matching process. The simplest way to 
do this is to imagine performing the matching slightly off-threshold, i.e. with the external 
particles having a small velocity in the center of mass frame. As mentioned earlier, this 
does not affect the value of the NRQED coefficients. Which NRQED diagram must be kept- 
is then decided by counting powers of a and v, with v ~ a kept in mind for bound-state 
applications. 

For example, for the present purposes of computing the 0(ma 5 ) hyperfine splitting, 
we show in §3 that we require both C 4 and C 5 to next-to-leading order, 0(o 2 ). This implies 
both of these couplings must be matched to QED with an accuracy up to order a 2 v°. The 
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NRQED diagram proportional to c:^ is of order a 2 v°. The loop diagram involving the 
exchange of one Coulomb photon contributes to order 3 a 2 /v and cancels a similar term 
in the QED one loop vertex correction. At threshold, the 1/v contribution becomes 4 a 
1/A infrared divergence which, cancels a similar term in the QED diagram. Because the 
coefficients of the effective theory describes only high-energy virtual effects in QED, it is 
a general result that matching always produces infrared finite values for them. 

Notice also that all other NRQED loop graphs which could appear on the right-hand- 
side necessarily involve additional powers of the electron or positron velocity, v, and so give 
contributions to C 4 which are smaller than 0(oi 2 ). For example, since the Feynman rule 
for the emission of a soft photon is proportional to ep/rn = ev, the loop graph obtained 
by dressing a four-fermi interaction with a soft-photon line rather than a Coulomb line, 
contributes a term of order au 2 cf to the scattering amplitude, representing a correction 
to C 4 which is of order a 4 in any bound-state application. In this way it may be seen that 
only the Coulomb-exchange NRQED loop need be kept in Fig. (3). 
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Figure (4) 

Diagrams which contribute to two-photon annihilation contributions to c^. 

We next turn to the QED box graphs describing s-channel electron-positron annihila¬ 
tion into two virtual photons, Figure (4). The matching appropriate for these graphs gives 
[ 8 ], [9]: 

C4 1 ' ) (2-7 aim) = 0, (2-y ann) = (% — 2 In 2 + i7^j. (7) 

This time charge-conjugation invariance forbids a contribution to the spin-triplet operator, 

c 4 . 

3 This diagram is proportional to a 2 J d 3 k/((p 2 — k 2 +ie)(k—p) 2 )~a 2 /p where p=mv is the external momentum. 

4 We regulate all such infrared divergences by including a photon mass, A, into photon propagators. 
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Notice, in this last equation, that the one-loop contribution to C5, has both a real 
and imaginary part. This imaginary part causes the low-energy hamiltonian not to be 
hermitian. The resulting loss of unitarity in the time evolution is just what is required to 
describe the depletion of electrons and positrons due to their mutual annihilation into real 
photons. Since annihilation is a high-energy effect, it appears in NRQED as an effective 
four-fermi operator. This imaginary part may be used to compute the decay rate for 
positronium bound states by calculating the imaginary part it implies for the bound state 
energy eigenvalue, E. The decay rate for bound state is then given by the familiar relation, 
T = -2 lm(E). 
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Figure (5) 

The one-loop t-channel matching diagrams. 

We now turn to the contribution to C 4 and C 5 due to the /.-channel two-photon QED 
exchange graphs, given in Figure (5). Here the diagram proportional to C is meant to 
represent the vertex, c^\ when the diagrams are evaluated in a spin 1 state, and 1 when 
the diagram is evaluated in a spin 0 state. It is straightforward to verify that only the 
given NRQED graphs can contribute to these coefficients to 0(o 2 ') when v ~ a. 

In order to determine and separately, we must include all of the spin-independent 
NRQED diagrams in the matching. This is actually more work than is required purely for 
the purposes of calculating the hyperfine splitting, since only the difference c A ' f — c 2 f enters 
into this quantity. Since the spin-independent graphs cancel in this difference, it suffices 
to just compute spin-dependent graphs if one is strictly interested only in the hyperfine 
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splitting. However, we present here the separate matching for both C4 and C5, since these 
two coefficients must be known separately for other applications, such as for the complete 
0(ma 5 ) shift of the positronium energy levels. 

Evaluating the left-hand side of Fig. (5) (the QED graphs) for a spin-triplet electron- 
positron configuration gives, after a straightforward calculation: 
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Recall A is the infrared-regulating photon mass. 

Using the NRQED Feynman rules to calculate the right-hand-side of Fig. (5) gives 
the following contributions. The only spin-dependent NRQED diagram, part ( d ) of Figure 
(5), gives: 
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where the overall factor of 2 in the first line takes into account the two possible ways in 
which the NRQED diagram can be drawn. 

Similarly, diagram (/) of Fig. (5) gives: 
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All the other NRQED diagrams can be calculated in a similar manner. The final 
result for the sum of the NRQED diagrams, evaluated in a spin-triplet state is: 


(NRQED) s= i - -2c { ^\t-ch) + a 


Solving for C 4 (f-ch) then gives the result: 
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An identical procedure applies to the S = 0 state. The QED graphs are then found 
to give: 

(QED) s =o = « 2 
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whereas the other NRQED diagrams are left unchanged because they are spin independent. 
The sum of the NRQED diagrams is then found to be 
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The complete one-loop contributions to the coefficients C 4 and C 5 are then given by 
combining the results from all three classes of graphs: 
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With the NRQED lagrangian now in hand, we next proceed with the determination 
of which graphs can contribute to the 0(ma 5 ) hyperfine splitting in positronium. 

3. Power Counting 

The essence of any effective field theory is its powercounting rules, since these are 
what permits the systematic calculation of observables to any order in small ratios of 
scales. Unfortunately, power counting in NRQED is slightly more complicated than in 
many effective field theories because of the appearance of ‘ultrasoft’ photons. Recall that 
the charged particles in a QED bound state typically have momenta p ~ ma and energy 
E ~ ma 2 . Photons can therefore be emitted with momenta equal to either of these scales. 
Photons having momenta, k, (and energy, uj) of order ma (which we will refer to as ‘soft’) 
do not pose any problems for powercounting, but those having k, u ~ ma 2 — the ‘ultrasoft’ 
ones — do. Physically, such ultrasoft photons represent the effects of retardation in the 
effective theory. 

Fortunately, these ultrasoft photons have wavelengths which are also large compared 
to the size of the bound state, and so their effects can be organized into slightly more 
complicated powercounting rules using what amounts to a multipole expansion in their 
couplings to the charged particles. In the final analysis, this multipole expansion introduces 
extra suppression by powers of a into interaction vertices involving ultrasoft photons and 
their contribution to the hyperfine splitting starts at order ma 6 [2], 

3.1) General Powercounting Rules 

When the dust settles, the powercounting result as applied to positronium has an 
appealingly simple form [2], Consider computing a contribution to a positronium bound- 
state observable using the NRQED lagrangian, Eq. (1), in old-fashioned nonrelativistic 
Rayleigh-Schrodinger perturbation theory. Each term in this perturbation series can be 
expressed graphically, with time moving to the right, say, and with internal lines denoting 
the particles which contribute to the sum over intermediate states. 
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For NRQED, we use Coulomb gauge, which is the most efficient gauge for the study 
of nonrelativistic systems. In this gauge, there are two types of photon propagator, one 
corresponding to Coulomb photons (which are represented by dotted lines in the Feyn¬ 
man diagrams) and one corresponding to transverse photons (represented by wavy lines). 
The Coulomb photons are soft whereas the transverse photons can be either soft or ultra- 
soft. Soft transverse photons are represented by slanted lines in old-fashioned perturbation 
theory whereas soft transverse photons are represented by vertical lines (see [2] for more 
details). As mentioned above, the ultra-soft contributions can be shown to be beyond the 
order we are interested in this work, and we will therefore consider only soft transverse 
photons. When only soft photons are present - i.e. when retardation effects are neglected 
- all NRQED diagrams reduce to a set of instantaneous interactions, which can be Fourier 
transformed to local interactions in coordinate space separated by electron-positron prop¬ 
agators. 

There are three quantities which determine the size of the contribution of any such 
graph to bound-state observables [2]. For any graph define N to be the number of electron- 
positron propagators separating the instantaneous interactions. 5 The other two quantities 
are related to the powers of a and 1/m which appear in each of the coupling constants, c % . 
of the NRQED lagrangian. Define, then, k and n to be, respectively, the total number of 
powers of 1/m (or v) and a which appear in the vertices of the graph of interest. 

Suppose N, k and n are known for any particular NRQED graph. Then the contri¬ 
bution of this graph to the energy-level shift in positronium depends on m and a (modulo 
logarithms of a) through the combination ma p , where [2]: 

p=l + K + n — N. (18) 

Although the quantity N enters this expression with a negative sign, inserting additional 
interactions (and so increasing N ) typically involves sufficient additional vertices to ensure 
that the net contribution to p increases and that the contribution from the diagram is 
therefore suppressed. The only exception to this statement is repeated insertions of the 

° N is denoted N T op in ref. [2]. 
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Coulomb interaction, which increase both N and n by one but leaves k unchanged (recall 
the Coulomb interaction does not contain any powers of inverse mass). Therefore, adding 
any number of Coulomb interactions to a given diagram leaves the value of p unchanged, 
indicating that one cannot perturb in the Coulomb interaction, which must be summed up 
to all orders. This is accomplished by using Schrodinger wavefunctions for the external lines 
of the bound state diagrams and the nonrelativistic Schrodinger-Coulomb propagator for 
intermediate states. On the other hand, it is easy to see that adding any other interaction 
increases the value of p since they always contain powers of 1/m. There is therefore only 
a finite number of graphs which can contribute for any positive choice for p. 

3.2) Applications to HyperGne Splitting 

We may now apply these powercounting arguments to the hyperfine splitting. An 
important simplifying feature appears in this specific application, since the hyperfine split¬ 
ting compares the energies of two states having different net spin. As a result it suffices 
to consider only graphs for which at least one of the vertices involves a spin-dependent 
coupling. 

To proceed, we start by recapping the powercounting which leads to the graphs which 
contribute to the 0(ma 4 ) hyperfine splitting [4], In this case we require p — 4, so k + n — 
N = 3. Consider first the case N — 0. In this case we require all possible graphs containing 
a single instantaneous photon exchange, subject to the following two conditions: (i) at least 
one vertex is spin dependent; and (ii) k + n = 3. The five classes of graphs which satisfy 
these two conditions are displayed in Figure (6). The coupling which appears in part (a) 
of this figure is ci, which appears twice, so inspection of Eq. (2) shows that the tree-level 
contribution, c^\ gives k = 2 and n — 1, as required. The same is true for parts (b) 
- which is linear in — and (c) — which is proportional to Cg 0 ^. Finally, parts (d) 
and (e) involve the four-fermi couplings, C4 and C5, the largest of which also arises at 
tree-level, proportional to a/m 2 , again giving n = 2, n = 1. It follows that these graphs 
contribute to 0(ma 4 ) provided we use the leading-order contribution to their couplings. 
It is noteworthy that only the graphs of parts (a), (d) and (e) of Figure ( 6 ) contribute to 
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Figure (6) 


The NRQED graphs which contribute to the hyperfine structure at order ma 4 . 

the hyperfine splitting of the ground state, since all of the others vanish when evaluated 
in an s-wave configuration. 

It now remains to consider the case N > 0. Having N = 1 implies including one more 
interaction in addition to the diagrams shown in Figure (3). It can be easily verified that 
any additional interaction (excluding, of course, the Coulomb interaction) increases the 
combination k + n by at least 3. For example, one can add a transverse photon coupled to 
two dipole vertices, introducing the square of the coefficient c'j 0 '* = qe/ (2m). This increases 
k by 2 and n by 1. The net change in p is therefore 3 — 1 = 2 and the diagram contributes 
only to order ma 6 . Similarly, adding a Coulomb photon connected to a Coulomb and 
Darwin vertex increases k by 2 and n by 1, whereas adding a relativistic kinetic vertex 
doesn’t change n but increases n by 3. Adding other interactions necessarily leads to, at 
best, 0(ma 6 ). 

This last argument has immediate implications for calculating the 0(ma 5 ) hyperfine 
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structure. To this order all of the relevant graphs must still have N = 0. The next- 
to-leading result is therefore obtained from the same graphs, Figure (6), but using the 
next-to-leading order — i.e. 0(a 2 ) — contributions to the coefficients, c t . We next show 
in more detail how these corrections are computed for the case of the coefficients of the 
four-fermi interactions, c^ and . 


4. Results 


To compute the hyperhne splitting we evaluate the graphs of Figure (6) using Coulomb 
wavefunctions to describe the initial and final electron-positron lines. Consider, first, the 
hyperhne splitting for s-wave states. In this case only graphs (a), (d) and (e) of Fig. (6) 
contribute since the other two graphs contain one vector, a, which can’t be dotted into 
any other vector if £ = 0. We find 
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(19) 


where n and £ — 0 are the principal and orbital quantum number of the positronium state 
of interest, and S — 0 or 1 is the net intrinsic spin of the e + e“ state. 

Using the formula T = —2 Im (E), we obtain from SE n (e) the 0(ma 5 ) decay rate of 
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the s-wave state of parapositronium (S — 0): 


r(n, a 5 ) 


ma 5 


( 20 ) 


In the rest of the paper, we concentrate on the hyperfine splitting so we drop the imaginary 
part contained in SE n (e). Adding the contributions of diagrams (a), (d) and (e) of Fig. (6), 
and taking the difference between S = 1 and S — 0 finally gives: 
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( 21 ) 


in agreement with standard results [6] [7]. 

The hyperfine splitting for arbitrary quantum numbers is computed by modifying the 
previous calculation in two ways. First, diagram (a) of Fig. (6) must be computed for £ ^ 0 
states. Notice that diagrams (d) and (e) of this figure are nonzero only for s-wave states 
since they represent contact interactions. Secondly, diagrams (6) and (c), which contribute 
only to £ 0 states, must be computed. Since the calculation of these latter contributions 

is presented elsewhere [4], we present here just the final result: 
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( 22 ) 


where j is the total angular momentum quantum number, including electron spin, and the 
coefficients Cji are given explicitly by: 


Cji = 


4^+5 

2(^+l)(2£+3) 

1 

-M+l 

2 £( 2 £—l) 


if j — £ + 1 
if j = i 
if j = £-l 


(23) 


To summarize, we see that the hyperfine splitting to O (mo 5 ) is hardly more difficult 
to obtain in NRQED than is the 0(ma A ) result. The only extra effort required is obtaining 
the complete matching of all spin-dependent effective operators to next-to-leading order in 
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a. We have performed the required 0(o 2 ) matchings for those four-fermi operators which 
had not been previously given in the literature. Furthermore, results for the hyperfine 
splitting for general n and £ are obtained with very little effort. 
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NRQED Feynman Rules. 
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Spin-Seagull Vertex 
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